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Abstract 

The method of "Isotropic Entanglement" (IE), inspired by Free Probability Theory and 
Random Matrix Theory predicts the eigenvalue distribution of quantum many-body (spin) systems 
with generic interactions. At the heart is a "Slider", which interpolates between two extrema by 
matching fourth moments. The first extreme treats the non-commuting terms classically and the 
second treats them isotropically. Isotropic means that the eigenvectors are in generic positions. 
We prove Matching Three Moments and Slider Theorems and further prove that the interpolation 
is universal, i.e., independent of the choice of local terms. Our examples show that IE provides an 
accurate picture well beyond what one expects from the first four moments alone. 
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I. ELUSIVE SPECTRA OF HAMILTONIANS 



Energy eigenvalue distributions or the density of states (DOS) are needed for calculating 
the partition function jl|, p. 14]. The DOS plays an important role in the theory of solids, 
where it is used to calculate various physical properties such the internal energy, the density 
of particles, specific heat capacity, and thermal conductivity {2!, Q . Quantum Many-Body 
Systems (QMBS) spectra have been elusive for two reasons: 1. The terms that represent 
the interactions are generally non-commuting. This is pronounced for systems with ran- 
dom interactions (e.g., quantum spin glasses [ig!, p. 320] 17, 18|). 2. Standard numerical 
diagonalization is limited by memory and computer speed. Calculation of the spectrum of 
interacting QMBS has been shown to be difficult jsl. 

An accurate description of tails of distributions are desirable for condensed matter physics. 
Though we understand much about the ground states of interacting QMBS js - 12|, eigenvalue 
distributions are less studied. Isotropic Entanglement (IE) provides a direct method for 
obtaining eigenvalue distributions of quantum spin systems with generic local interactions 



and does remarkably wel 
entanglement generation 



in capturing the tails. Indeed interaction is the very source of 
Section 2.4.1] 14] which makes QMBS a resource for quantum 
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computation 15| but their study a formidable task on a classical computer. 

Suppose we are interested in the eigenvalue distribution of a sum of Hermitian matrices 
M = YliLi^i- general, Mj cannot be simultaneously diagonalized, consequently the 
spectrum of the sum is not the sum of the spectra. Summation mixes the entries in a very 
complicated manner that depends on eigenvectors. Nevertheless, it seems possible that a 
one-parameter approximation might suffice. 

Though we are not restricted to one dimensional chains, for sake of concreteness, we 
investigate interacting (i-dimensional quantum spins (qudits) on a line with generic inter- 
actions. The Hamiltonian is 

N-l 

H=Y^ ^d^-' ® ^'r- ,l+L~l ® Id^-i-(L-i) , (1) 
1=1 

where the local terms Hi ... i^l_i are finite x d^ random matrices. We take the case of 
nearest neighbors interactions, L = 2, unless otherwise specified. 

The eigenvalue distribution of any commuting subset of H such as the terms with / odd 
(the "odds") or / even (the "evens") can be obtained using local diagonalization. However, 
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Figure 1: Odd and even summands can separately be locally diagonalized, but not the sum. The 
overlap of the two subsets at every site generally requires a global diagonalization. 

the difficulty in approximating the full spectrum of H = -f^odd + -^even summing the 
odds and the evens because of their overlap at every site. 

The intuition behind IE is that terms with an overlap, such as Hi^i+i and -f^/+i,/+2) intro- 
duce randomness and mixing through sharing of a site. Namely, the process of entanglement 
generation introduces an isotropicity between the eigenvectors of evens and odds that can 
be harnessed to capture the spectrum. 

Random Matrix Theory (RMT) often takes advantage of eigenvectors with Haar measure, 
the uniform measure on the orthogonal/unitary group. However, the eigenvectors of QMBS 
have a more special structure (see Eq. [TT]) . 

Therefore we created a hybrid theory, where we used a finite version of Free Probability 
Theory (FPT) and Classical Probability Theory to capture the eigenvalue distribution of 
Eq. [U Though such problems can be QMA-complete, our examples show that IE provides 
an accurate picture well beyond what one expects from the first four moments alone. The 
Slider (bottom of Figure [2]) displays the proposed mixture p. 

II. THE METHOD OF ISOTROPIC ENTANGLEMENT 
A. Overview 

We propose a method to compute the "density of states" (DOS) or "eigenvalue density" 
of quantum spin systems with generic local interactions. More generally one wishes to 
compute the DOS of the sum of non-commuting random matrices from their, individually 
known, DOS's. 

We begin with an example in Figure |3l where we compare exact diagonalization against 
two approximations: 

• Dashed grey curve: classical approximation. Notice that it overshoots to the right. 
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add random variables = 
Classical convolution of densities 



add random variables = 
Classical convolution of densities 



Section lll-B,C 



Isotropic (or free) convolution = 
dv"° = dv„ ffl, dv, 

A ISO B 



Several options for adding 
random variables 



Section IV 



Isotropically Entangled 
(1-p) dv'=° + p dv^ 




Isotropic 
p = 



Section lll-D 



Section lll-A 



Classical convolution = 
dv"^ = dv„ * dv„ 



Quantum convolution = dv'' = dv. ffl„dv„ = 

A q D 

Isotropically Entangled 



Classical 
p = l 



Figure 2: The method of Isotropic Entanglement: Quantum spectra as a convex combination of 
isotropic and classical distributions. The Slider (bottom) indicates the p that matches the quantum 
kurtosis as a function of classical (p = 1) and isotropic (p = 0) kurtoses. To simplify we drop the 
tensor products (Eq. [5]) in the local terms (ellipses on top). Note that isotropic and quantum 
convolution depend on multivariate densities for the eigenvalues. 

• Solid grey curve: isotropic approximation. Notice that it overshoots to the left. 

• Solid black curve: isotropic entanglement (IE). 

• Dots: exact diagonalization of the quantum problem given in Eq. [H 

The classical approximation ignores eigenvector structure by summing random eigenvalues 
uniformly from non-commuting matrices. The dashed curve is the convolution of the prob- 
ability densities of the eigenvalues of each matrix. 
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Figure 3: The exact diagonalization in dots and IE compared to the two approximations. The title 
parameters are explained in the section on numerical results. 

The isotropic approximation assumes that the eigenvectors are in "general position" ; that 
is, we add the two matrices with correct eigenvalue densities but choose the eigenvectors 
from Haar measure. As the matrix size goes to infinity, the resulting distribution is the free 
convolution of the individual distributions [22]. 

The exact diagonalization given by red dots, the dashed and solid grey curves have exactly 
the same first three moments, but differing fourth moments. 

Isotropic Entanglement (IE) is a linear combination of the two approximations that is 
obtained by matching the fourth moments. We show that 1) the fit is better than what might 
be expected by four moments alone, 2) the combination is always convex for the problems 
of interest, given by < p < 1 and 3) this convex combination is universal depending on 
the parameter counts of the problem but not the eigenvalue densities of the local terms. 

Parameter counts: exponential, polynomial and zero. Because of the locality of generic 
interactions, the complete set of eigenstates has parameter count equal to a polynomial in 
the number of spins, though the dimensionality is exponential. The classical and isotropic 
approximations have zero and exponentially many random parameters respectively. This 
suggests that the problem of interest somehow lies in between the two approximations. 

Our work supports a very general principle that one can obtain an accurate representa- 
tion of inherently exponential problems by approximating them with less complexity. This 
realization is at the heart of other recent developments in QMBS research such as Matrix 
Product States 0,0], and Density Matrix Renormalization Group ^], where the state (usu- 



6 



Inputs 


^ - ^odd + ^even ^ ^^d, (3 


— > 


Output 


p, dv^^ 



Figure 4: Inputs and outputs of the IE theory. See section Hill for the definition of dv^^ . 

ally the ground state of ID chains) can be adequately represented by a Matrix Product 
State (MPS) ansatz whose parameters grow linearly with the number of quantum particles. 
Future work includes explicit treatment of fermionic systems and numerical exploration of 
higher dimensional systems. 



B. Inputs and Outputs of the Theory 

In general we consider Hamiltonians H = H^^^ "'"-^even' ^^lere the local terms that add 
up to -f^Qf^j (or -f^gven) fo^'^ ^ commuting subset. All the physically relevant quantities 
such as the lattice structure, A^, dimension of the spin d and the rank r are encoded in the 
eigenvalue densities. The output of the theory is a < p < 1 by which the IE distribution 
is obtained and du^^ serves as an approximation to the spectral measure. The inputs can 
succinctly be expressed in terms of the dimension of the quantum spins, and the nature of 
the lattice (Figure Hj). 



C. More Than Four Moments of Accuracy? 

Alternatives to IE worth considering are 1) Pearson and 2) Gram-Charlier moment fits. 
We illustrate in Figure [5] how the IE fit is better than expected when matching four 
moments. We used the first four moments to approximate the density using the Pearson 



fit as implemented in MATLAB and also the well-known Gram-Charlier fit [2J]. In [23 1 
it was demonstrated that the statistical mechanics methods for obtaining the DOS, when 
applied to a finite dimensional vector space, lead to a Gaussian distribution in the lowest 
order. Further, they discovered that successive approximations lead naturally to the Gram- 



Charlier series 2J] . Comparing these against the accuracy of IE leads us to view IE as more 
than a moment matching methodology. 

The departure theorem (Section IIVBI) shows that in any of the higher moments (> 4) 
there are many terms in the quantum case that match IE exactly. Further, we conjecture 
that the effect of the remaining terms are generally less significant. 
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4-moment matching methods vs. IE: binomial local terms 



4-moment matching methods vs. IE: Wishart local terms 
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Figure 5: IE vs. Pearson and Gram-Charlier 
III. SPECTRA SUMS IN TERMS OF PROBABILITY THEORY 

The density of eigenvalues may be thought of as a histogram. Formally for an m x m 
matrix M the eigenvalue distribution is 19|, p. 4] 20|, p. 101] 



^ m 
i=l 



(2) 



For a random matrix, there is the expected eigenvalue distribution [2l|, [22|, p. 362] 

1 



-E 



m 



Y,5{x-\ (M)) 



i=l 



(3) 



which is typically a smooth curve [l9i, p. 101]|20|, p. 115]. 

The eigenvalue distributions above are measures on one variable. We will also need the 
multivariate measure on the m eigenvalues of M: 

dvM {x) = The symmetrized joint density of the eigenvalues. 

Given the densities for M and M', the question arises: What kind of summation of 
densities might represent the density for M + M'? This question is unanswerable without 
further information. 

One might try to answer this using various assumptions based on probability theory. The 
first assumption is the familiar "classical" probability theory where the distribution of the 
sum is obtained by convolution of the density of summands. Another assumption is the 
modern "free" probability theory; we introduce a finite version to obtain the "isotropic" 
theory. Our target problem of interest, the "quantum" problem, we will demonstrate, prac- 
tically falls nicely in between the two. The "Slider" quantifies to what extent the quantum 
problem falls in between (Figure [2] bottom). 
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A. Classical 



Consider random diagonal matrices A and B of size m, the only randomness is in a 
uniform choice among the ml possible orders. Then there is no difference between the 
density of eigenvalue sums and the familiar convolution of densities of random variables, 



(4) 



Comment: From this point of view, the diagonal elements of A, say, are identically 
distributed random variables that need not be independent. Consider Wishart matrices 
where there are dependencies among the eigenvalues. To be precise let a G be 
a symmetric random variable, i.e.. Pa has the same distribution as a for all permutation 



matrices P. We write, A 



ai 



\ 



\ 



J 



= diag(a). Similarly for B. 



Comment: The classical convolution appears in Figure [2] in two different ways. Firstly, in 
the definition of A (or B) , the eigenvalues of the odd (or even) terms are added classically. 
Secondly, A and B are added classically to form one end of the Slider. 



B. Free and Isotropic 

Free probability j22i, is recommended] provides a new natural mathematical "sum" of 
random variables. This sum is computed "free convolution" denoted 

du^ = duA ffl duB- (5) 

Here we assume the random matrices A and B, representing the eigenvalues, have densities 
duA and duB- In the large m limit, we can compute the DOS of A + Q^BQ, where Q is a 
/3— Haar distributed matrix (see Table [T]). 

Comment: In this paper we will not explore the free approach strictly other than observing 
that it is the infinite limit of the isotropic approach (i.e., t — )■ cxd in Eq. E]). This infinite 
limit is independent of the choice of /3 (see Table [T]). 
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Table I: Various Haar matrices. 



We define an isotropic convolution. The isotropic sum depends on a copying parameter 
t and /3 (Table [T]). The new Hamiltonian is the isotropic Hamiltonian ( "iso"): 

H,,o = {A' ® li) + {\ ® B') Qp, (6) 

where Q/s is a /3— Haar distributed matrix, A = A' ®lt and B = lt® B' . For the copying 
parameter t = d, Hiso has the same dimension as H in Eq. [U however, t > d allows us 
to formally treat problems of growing size. We can recover the free convolution by taking 
the limit: limj_j.oo c^i^***"'-^'*-' = du^ . The effect of Qp is to spin the eigenvectors oilt® B to 
point isotropically with respect to the eigenvectors of A. We denote the isotropic eigenvalue 
distribution by 

omitting t and /3 when it is clear from the context. 

Comment: In Eq. |6l the If and i? in Ij (g) B, can appear in any order. We chose this 
presentation in anticipation of the quantum problem. 

Comment: In this paper we primarily consider t to match the dimension of H. 

C. Quantum 

Let du'^ denote the eigenvalue distribution for the Hamiltonian in Eq. [H This is the 
distribution that we will approximate by du^^. In connection to Figured] the Hamiltonian 
can be written as 

^ = ^odd + ^even = I ® ^M+i ® ^ + ^ ® ^M+i ® I- (») 

l=l,3,5,--- l=2,4,e,--- 
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We proceed to define a "quantum convolution" on the distributions duA and dOs, which 
is /3- dependent 



dui^P^ = dp A ffl« dpB- 



(9) 



In general, without any connection to a Hamiltonian, let dvA and dvB be symmetric 
measures on W^'^ . We define dz/'^*^^) to be the eigenvalue distribution of 



H = A + Q-'BQ 



11 



where Q„ 



Qf ^ with 



(10) 



(iV-l)/2^(0) 



h and Qf ) = I, 



AT odd 



N even 



(11) 



and each Qj-*"* is a /3— Haar measure orthogonal matrix of size d"^ and 1^ is an identity 
matrix of size d. 



Comment: A, B and Qq are d^ x d^ . 
Comment: In our examples given in this paper, we assume the local terms are independent 
and identically distributed (iid) random matrices, each of which has eigenvectors distributed 
with /3— Haar measure. 

The tensor product in (fTTj) succinctly summarizes the departure of the quantum case from 
a generic matrix as well as from the classical case. First of all the number of parameters in 
Qq grows linearly with whereas in Q it grows exponentially with A^. Second, the quantum 
case possesses isotropicity that makes it different from the classical, whose eigenvectors are 
a point on the orthogonal group (i.e., the identity matrix). 

Comment: General /3's can be treated formally In particular, for quantum mechan- 
ical problems /3 is taken to be 1 or 2 corresponding to real and complex entries in the local 
terms. /3 = 4 corresponds to quaternions. 
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Definition. The Hadamard product of two matrices Mi and M2 of the same size, denoted 
by Ml o M2, is the product of the corresponding elements. 

Lemma 1. The elements of Qq defined in Eq. IIJ\ are (dependent) random variables with 
mean zero and variance . 

Proof. Here expectations are taken with respect to the random matrix Q„ which is built 



from local Haar measure matrices by Eq. [TTJ The fact that E ( Q 



E ( gf ) 



follows from the Haar distribution of local terms. Thus E {Qq) = by independence of 
Q^^"* and Qg^^- Further, each element in Qq involves a dot product between columns of 
gS^) and gf ^ In every given column of Qq any nonzero entry is a distinct product of 
entries of local Q's (see Eq JTT]) . For example the expectation value of the 1,1 entry is 
IE [qI^^ <l)j^ (ii^i Q^^i^ = E {(li^<l^j^^ E Because of the Haar measure of the local 

terms, this expectation is zero unless i = j. We then have that 



E (g, o g,) = E (gj^) o g(^)) e (gf ^ o g 

= d'^JaN, 

2 



1 



Ir 



AT odd 
N even 



(12) 



where Ji = i x i matrix of all ones. We use facts such as (J-i/i) = (Ji/i), (Ji/i) 
(Jj2/i^) and the variance of the elements of an i x -i /3— Haar matrix is 1/i. 



® {Ji/i) 



□ 



IV. THEORY OF ISOTROPIC ENTANGLEMENT 



A. Isotropic Entanglement as the Combination of Classical and Isotropic 



We create a "Slider" based on the fourth moment. The moment of a random variable 
defined in terms of its density is = J x^dvM- For the eigenvalues of an m x m random 
matrix, this is — ETrM*^. In general, the moments of the classical sum and the free sum 
are different, but the first three moments, mi, m2, and 7713 are theoretically equal 22|, p. 
191]. Further, to anticipate our target problem, the first three moments of the quantum 
eigenvalues are also equal to that of the iso and the classical (The Departure and the Three 
Moments Matching theorems in Section [iVBp . These moments are usually encoded as the 
mean, variance, and skewness. 
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We propose to use the fourth moment (or the excess kurtosis) to choose a correct p from 
a shding hybrid sum: 



du'^ ^ dv^^ = pdv^ + (1 - p)dv'''' (13) 

Therefore, we find p that expresses the kurtosis of the quantum sum (7I) in terms of the 
kurtoses of the classical (7I) and isotropic (72"*") sums: 

il = inl + {i-p)ir^ P=^$^o- (14) 

12 12 

Recall that the kurtosis 72 = where is the variance. Hence kurtosis is the correct 
statistical quantity that encodes the fourth moments: 



ml = — ETr {A + U^BU)\ mf = ^ETr {A + Q^BQ)\ ml = ^ETr {A + Q'^BQ^Y , 

(15) 

where 11 is a random uniformly distributed permutation matrix, Q is a /3— Haar measure 
orthogonal matrix of size d^ , and Qg is given by Eq. [Ill Unless stated otherwise, in 
the following the expectation values are taken with respect to random eigenvalues A and B 
and eigenvectors. The expectation values over the eigenvectors are taken with respect to 
random permutation U, fi—Haar Q or Qq matrices for classical, isotropic or quantum cases 
respectively. 



B. The Departure and The Matching Three Moments Theorems 

In general we have the i^^ moments: 



^SO ^ l^r^^ (A + Q^BQX 
m ^ ' 



q 

m; 



-KIiiA + Q^BQ,)\ and 



m 

jN 



ml = — ETr (A + n^5n)' 



where m = d . If we expand the moments above we find some terms can be put 
in the form ETr (^A"^'^ B"^^ Q and the remaining terms can be put in the form 
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ETr |. . . Q^B-^Q,A-^Q^B-^Q, . . .} . The former terms we denote non- departing; the re- 
maining terms we denote departing. 
For example, when i = 4, 



m 

rni = — E i Tr 



m 



^E |Tr ^A" + AA^Q^BQ + AA^Q^B^Q + AAQ^B^Q + 2 (AQ^BQ)^ + B^] } (16) 
+ AA^Q'^BQg + AA^Q^B^Q, + AAQ^B'^Qg + 2 (AQjBQq)^ + B 



rrf. 



— E (Xr + AA^lfBIi + AA^lf B'^Ii + AAlf B^Ii + 2 (AH^Bn)^ + 5^1 1 
m I L — — J J 



where the only departing terms and the corresponding classical term are shown as underlined 
and bold faced. 

Theorem. (The Departure Theorem) The moments of the quantum, isotropic and clas- 
sical sums differ only in the departing terms: ETrj. . . Q^B-^Q,A-^Q^B-^Q, . . .} . 

Proof. Below the repeated indices are summed over. If A and B are any diagonal matrices, 
and Q, is Q or Qq or 11 of size m x m then E (g|^) = 1/m , by symmetry and by Lemma 
[H respectively. Since the ETr (AQ^-BQ,) = E(g|^aj6j), where expectation is taken over 
randomly ordered eigenvalues and eigenvectors; the expected value is (^) E(aj6j) for 
any t or j. Hence, ^ETr (AQ^Sg.) = E{a^bj) = E{ai)E{bj), which is equal to the 
classical value. The first equality is implied by permutation invariance of entries in A and 
B and the second equality follows from the independence of A and B. □ 



Therefore, the three cases differ only in the terms ^ETr (AQ^BQ) , ^ETr (AQ^BQq) 
and ^ETr {AU^BUY in Eq. M 

Theorem. (The Matching Three Moments Theorem) The first three moments of the 
quantum, iso and classical sums are equal. 

Proof. The first three moments are 

mj'^ = ^ETr {A + B) 
mi'^ = i^ETr {A + QlBQ.) ' = ^ETr {A' + 2AQ^BQ. + B') (17) 



m 



(•) 



^ETr {A + Q^BQ.f = ^ETr {A' + SA^Q^BQ. + SAQ^B^Q. + B^) 



where Q, is Q and Qq for the iso and the quantum sums respectively and we used the 
familiar trace property Tr(MiM2) = Tr(M2Mi). The equality of the first three moments of 
the iso and quantum with the classical follows from The Departure Theorem. □ 
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number of 
sites 


number of odds 
or evens {N odd) 


number of odds 
or evens {N even) 


size of 
H 


Number of 
copies 




n — 2 


k - ^ k - 
odd 2 ' '^even 2 


m = 


t 



dimension of qudits 


size of local terms 


jth 22ioment 


Jth cumulant 


mean 


variance 


skewness 


kurtosis 


d 


n = d'^ 


mi 


1^1 






71 


72 



Table II: Notation 



Furthermore, in the expansion of any of the moments > 4 all the non-departing terms are 
exactly captured by IE. These terms are equal to the corresponding terms in the classical and 
the isotropic and therefore equal to any linear combination of them. The departing terms in 
higher moments (i.e.,> 4) that are approximated by IE, we conjecture are of little relevance. 
For example, the fifth moment has only two terms (shown in bold) in its expansion that are 
departing: 



m5 



i^ETr {a!> + 5A^Q^BQ. + SA^Q^B^Q. + SA^Q^B^Q. + 5A (AQ^BQ.)' + 



5 (AQ^BQ.)^ Q^BQ. + SAQ^B^Q, + fi^ 



By the Departure Theorem the numerator in Eq. [T3] becomes. 



(18) 



ll - 72 = - 

m 

and the denominator in Eq. [14] becomes. 



E<^Tr 



{AQ];BQ,f-{AQ^BQy 



72 - iT 



2 E |Tr [(An^^n)' - {aq^bqY'^ } 



(19) 



(20) 



where as before, Q is a /3— Haar measure orthogonal matrix of size m = , Q 



Qi^^ j Qf' given by Eq. [TT]and denote the fourth cumulants. Therefore, evaluation 
of p reduces to the evaluation of the right hand sides of Eqs. [19] and [20l 

Below we do not want to restrict ourselves to only chains with odd number of sites and 
we need to take into account the multiplicity of the eigenvalues as a result of taking the 
tensor product with identity. It is convenient to denote the size of the matrices involved by 
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m = = tn^^ where n = d? and k = ^^^^ and t is the number of copies (Section IIII Bl and 
Table [HI). 

C. Distribution of A and B 

The goal of this section is to express the moments of the entries of A and B (e.g., 
and m^]) in terms of the moments of the local terms (e.g for odd local terms m!^'^'^^ m'^^'^). 
Note that A and B are independent. The odd summands that make up A all commute and 
therefore can be locally diagonalized to give the diagonal matrix A (similarly for B), 

N~2 

A = ^ I®Ai(g)I (21) 

j=l,3,-- 

B = ^ I (g) Aj (g) I, 

j=2,4,--- 

where Aj are of size d'^ and are the diagonal matrices of the local eigenvalues. 

The diagonal matrices A and B are formed by a direct sum of the local eigenvalues of odds 
and evens respectively. For open boundary conditions (OBC) each entry has a multiplicity 
given by Table lllll 



OBC 


N odd 


N even 


A 


d 


1 


B 


d 


d^ 



Table III: The multiplicity of terms in A and B for OBC. For closed boundary conditions there is 
no repetition. 

Comment: We emphasize that A and B are independent of the eigenvector structures. In 
particular, A and B are the same among the three cases of isotropic, quantum and classical. 

We calculate the moments of A and B. Let us treat the second moment of A {B is done 
the same way). By the permutation invariance of entries in A 
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m^^E [al) = E {x^ + ■ ■ ■ + A^f 

= E[k{X^) +A;(A;- 1)AWa(')] (22) 

where expectation is taken over randomly chosen local eigenvalues, m^^^ = E (Af ) and 
^odd ^ ]g (A^Aj) for some uniformly chosen i and j with z ^ j. The permutation invariance 
assumption implies E (a^) = E {af) for alH = 1 ■ ■ ■ m. 

Comment: The key to this argument giving is that the indices are not sensitive to 
the copying that results from the tensor product with 1^ at the boundaries. 

Next we calculate the correlation between two diagonal terms, namely m^^ = E (aiaj) 
for i j. We need to incorporate the multiphcity, denoted by t, due to the tensor product 
with an identity matrix at the end of the chain. 



' E A<;' + ... + Af - j: (A!;> + ... + A!f) W23) 



mfi = — -E 

m{m — 1) 



= k{k -1)E {Xy + k {prob (A^) E (A^) + prob (A1A2) E (A1A2)} 

= k{k-l) mf"^ + — ^ I {tn'-' - 1) mOdd + (^^^-1 _ 1)) ^oddj 

where, prob (A^) = ^"^^"-"^ and prob (A1A2) = m-i'^^ ■ Similarly for B. 



D. Evaluation and Universality of p = l_ 



Recall the definition of p; from Eqs. [HI [19] and [20] we have. 



ETr (AU^BU) - ETr (AQ^BQ,) 

l-p= ^ S ^ ' . (24 

ETr {AW^ BUf -ETi{AQ'^BQ f ^ ' 



The classical case 



1 1 

^ETr (An^^n)' = -E ^ a\h\ = E (a^) E (fe^) = m^mf . (25) 



m m 

i=l 



Comment: Strictly speaking after the first equality we must have used Bt^. instead of 6^ 
but we simplified the notation as they are the same in an expectation sense. 
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moments 



m2mii 



(mil)' 



expectation values 



13+2 



m(m/3+2) 



E ( ki.igi,: 



m(m/3+2) 



E (91,1^1,292,192,2) 



/3 

m (m/3+2)(m-l) 



/3(n-l)+2 
n(n/3+2)(n-l) 



count 



2m^ (m — 1) 



w? (m — 1)' 



Table IV: The expectation values and counts of colliding terms in Q when it is /3— Haar distributed. 
In this section we use the first row; we include the last three rows for the calculations in the 
appendix. 

The general form for the denominator of Eq. [21] is 



-ETr 



m 



1 



{Ali'BIi) - {AQ^BQ) = -E {a^b^ - aiakb.bp {qjiq.kqpkqpi)] ■ (26) 



m 



It's worth noting that the arguments leading to Eq. |30]hold even if one fixes A and B and 
takes expectation values over 11 and a permutation invariant Q whose entries have the same 
expectation value. The right hand side of Eq. [30] is a homogeneous polynomial of order two 
in the entries of A and B] consequently it necessarily has the form 



-ETr 



m 



{AlfBYi) ' - {AQ^BQ) ' =ci{B,Q)m^ + C2 {B, Q) m^^ 



but Eq. [26] must be zero for A 
— C2, allowing us to factor out (? 
invariance of B implies, 



/, for which 



m 



1,1 



1. This implies that Ci 



i^. Similarly, the homogeneity and permutation 



-ETr 



[AIl^BIl) ' - [AQ^BQ) ' = (m^ - 1) [D^ (Q) mf + (Q) J 



The right hand side should be zero for B = I, whereby we can factor out ( 



rrio — m 



1,1) 



-ETr 



m 



{AU^BU)' - {AQ^BQY = {m^ - mf^,) (mf - mf J / (Q) 



(27) 



where = E (af ), mf = E , and m^^ = E (a^, aj) , mf ^ = E (6j, bj). Moreover / (Q) is 
a homogeneous function of order four in the entries of Q. To evaluate / (Q), it suffices to let 

th 

A and B be projectors of rank one where A would have only one nonzero entry on the i 
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position on its diagonal and B only one nonzero entry on the j position on its diagonal. 
Further take those nonzero entries to be ones, giving m^^^ = mfi = and = mf = 1/m, 



— ETr 

m 



{AlfBnY-{AQ^BQY]=^J{Q) 



But the left hand side is 



(28) 



-ETr 



m 



{AU^BuY-{AQ^BQy 



1 

m 

1 I 1 



m m 



1 f 1 



where, we used the homogeneity of Q. Consequently, by equating this to / {Q) jw? ^ we get 
the desired quantity 



/(g) = {l-mE (gj)} 



Our final result Eq. [27] now reads 



-ETr 



— m 



(29) 



The same calculation where each of the terms is obtained separately yields the same result 
(Appendix). In this paper p is formed by taking Q to have a /3— Haar measure. Expectation 
values of the entries of Q are listed in the Table HVl 

We wish to express everything in terms of the local terms; using Eqs. |22] and |23] as well 

as tn'' = m, 



mn — m 



1,1 



tk (n — l)n'' ^ 

m — 1 
tk {n — 1) n^~'' 

m — 1 



odd ™odd 



1,1 



[m. 



even ^even\ 



1,1 ; ' 



giving 
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-E 



m 



Tr {AU^BUY - Tr {AQ^BQY 



odd 



fcm (n — 1 



i^) {l-mE(4)}. 



(30) 



We now proceed to the quantum case where 



we 



need to evaluate 



-E 



{AU^BUY -TT{AQ^BQ^y 



In this case, we cannot directly use the tech- 
niques that we used to get Eq. [30] because Qq is not permutation invariant despite local 
eigenvectors being so. Before proceeding further we like to prove a useful lemma (Lemma 
E]). Let us simplify the notation and denote the local terms that are drawn randomly from 
a known distribution by -ffz,z+i = H^''^ whose eigenvalues are A; as discussed above. 

Recall that A represents the sum of all the odds and Q^^BQq the sum of all the evens, 



^odd = XI ^ ® ^^'^ ® ^' ^^^^ ^even = J] ^ ® ^^'^ ® ^' 

«=1,3,5,-- Z=2,4,6,--- 



Hence, the expansion of ^E |^Tr (^AQ^ BQq^ amounts to picking an odd term, an even 
term, then another odd term and another even term, multiplying them together and taking 
the expectation value of the trace of the product (Figure [6]). Therefore, each term in the 
expansion can have four, three or two different local terms, whose expectation values along 
with their counts are needed. These expectations are taken with respect to the local terms 
(dense (P x d"^ random matrices). 

The expectation values depend on the type of random matrix distribution from which the 
local terms are drawn. The counting however, depends on the configuration of the lattice 
only. We show the counting of the number of terms, taking care of the boundary terms for 
an open chain, along with the type of expectation values by which they need to be weighted: 

For odd {k odd terms and k even terms) 

Four i/W's -.k^k- if rf^-^^ETr , mi G {5, ■ ■ ■ , 8} 

Three i/W's : 2k^ {k - 1) =^ (/^-"^ETr ([/^^'^j^) ETr U2 G {4,5,6} 

Two ifW's: {k - if Not Entangled d^-^ |ETr ([i^^'^]^) }^ 
Two i/(-)'s : {2k - 1) Entangled ^ d^~^W£i [(ifC) l) (l O ® l) (l ® 

(31) 
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Q'BQ 

A 
Q'BQ 



can be visualized as picking an 



Figure 6: The terms in the expansion of Tr (^AQ^BQ^^ 
element from each row from top to bottom and multiplying. Each row has k of the local terms 
corresponding to a chain with odd number of terms. Among k"^ terms roughly A;^ of them differ 
among the classical, isotropic and quantum cases (See Eqs. [511 and [52]) . An example of such a 
choice is shown by diamonds. 

For even {k odd terms and k — 1 even terms) 

Four //W's -.kik- if {k - 2) ^ rf^-^^ETr (ifC))^ , mi G {5, ■ ■ • , 8} 
Three //W's : k {k - 1) {2k - 3) =^ ^^-"^ETr ([i^^'^]^) ETr {H^^^^f , U2 G {4,5,6} 
Two //W's : (A; - 1) (A; - 2) Not Entangled ^ d^-^ |ETr ([H^^^f^ }^ 

Two iJ(-)'s : 2 (A; - 1) Entangled ^ rf^-^ETr [{H^^^ ® l) (l O (iJ(') ® l) (l ® H'^^+^^)] 

(32) 

Here ui and U2 indicate the number of sites that the local terms act on (i.e., occupy). 



Therefore, — E 



TT{AQlBQ^y 



is obtained by multiplying each type of terms, weighted by 



the counts and summing. For example for ui = 5 and U2 = 3, when is odd. 



iE 



Tr [AQ^BQ.Y] = ^ [d^-^e {k - ifKIv (if«) ^ 



2A;2 {k - 1) rf^-^ETr ( [ifW] ETr (if«) ' + (A: - 1)' d^-^ jffiTr ( ') } + (33) 
(2A; - 1) c/^-3ETr [(i7« ® l) (l if('+^)) (i7« ® l) (l ® if('+i))] } 

and similarly for N even. 
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i^E [Tr {AQ^^BQ.Y] = (k - 1) (k - 2) d^'^ETi + 

k {2k - 3) rf^-^ETr ( [ifW] 2 j ETr ' + (A; - 2) d^-^ |ETr ( }' + (34) 
2rf^-3ETr ® I) (I ® ® I) (I ® } . 

The expectation values depend on the type of random matrix distribution from which 
the local terms are drawn. We will give explicit examples in the following sections. In the 
following lemma, we use 

Lemma 2. In calculating the E,Tr(^AQ^BQqY if at least one of the odds (evens) commutes 
with one of the evens (odds) then the expectation value is the same as the classical expectation 
value. Further if the local terms have permutation invariance of eigenvalues then the only 
quantum expectation value that differs from classical is of Type II (see the proof and the 
diamonds in figure 

Proof. This can be shown using the trace property Tr (MP) = Tr (PM). In calculat- 
ing ETr (^H°'^'^Hp'"^"'H°'^'^H^"'^"'^ ; if any of the odd (even) terms commutes with any of the 
even (odd) terms to its left or right then they can be swapped. For example one gets 

ETr (//orf^i^e^^en^et^en^oddj ^ ^rj.^ (^jjodd jjodd jjeven jjeven^ ^J^j^j^ j^g^ ^J^g claSSical Value. 

Hence the only types of expectations that we need to worry about are 

H'^^) and 

Type I Type II 

now we show that with permutation invariance of the local eigenvalues the first type are 
also classical leaving us with the "diamond terms" alone (Fig. [6]). Consider a Type I term, 
which involves three independent local terms. 



-ETr 



[(n^2 &^M-a) (l(E)H(2') (xDl^jv_3 ) (1^2 ®I^jv-4 ) (l^S ^^(4) ®I^jv_5 )] 
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This follows immediately from the independence of H^^^ , which allows us to take its expec- 
tation value separately giving a /i and leaving us with 



^ETr 

m 



(Irf2 ® if ^ (I ® if ^ 



Therefore the only relevant terms, shown by diamonds in Fig. El are of Type II. As an 
example of such terms consider (here on repeated indices are summed over) 

d3 y «l«2,ilJ2 liP2,jlk2 J y 321'iMk'j, %2t'i,P2k'j, J J ' 

where the indices with subscript 2 prevent us from treating the two expectation values 
independently: H^^'> and ii*^^) overlap at the second site. The number of such terms is 
2k - 1, where k = □ 

Therefore, we have found a further reduction of the terms from the departure theorem, 
that distinguishes the quantum problem from the other two. Luckily and interestingly the 
kurtosis of the quantum case lies in between the classical and the iso. We emphasize that 
the only inputs to the theory are the geometry of the lattice (e.g., the number of summands 
and the inter- connectivity of the local terms) and the moments that characterizes the type 
of the local terms. 

Comment: The most general treatment would consider Type I terms as well, i.e., there is 
no assumption of permutation invariance of the eigenvalues of the local terms. This allows 
one to treat all types of local terms. Here we are confining to random local interactions, 
where the local eigenvectors are generic or the eigenvalues locally are permutation invariant 
in the expectation value sense. 

The goal is to find p by matching fourth moments 

ETr [AW^BuY - ETr {AQ^BQ^f 
ETr {AlFBIif - ETr {AQ^BQf 
for which we calculated the denominator resulting in Eq. [301 where E (lo'jjl'^) = ^^^"^^^j 
/3— Haar Q (Table [TV]) . If the numerator allows a factorization of the moments of the local 
terms as in Eq. [301 then the value of p will be independent of the covariance matrix (i.e., 
eigenvalues of the local terms). 
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Lemma. (Universality) p t— )■ p{N,d, P), namely, it is independent of the distribution of 
the local terms. 

Proof. We use a similar techniques as we did in the isotropic case. The general form for the 
numerator of Eq. [21] is (denoting Lemma [2] by L3) 



— ETr 

m 



rf3 

{2k - 1) 



- [{Qi'^iQi ® h) {h ® gr+Wig^+i)]'} (se) 

where the expectation on the right hand side is taken with respect to the local terms H^^^ 
and . The right hand side is a homogeneous polynomial of order two in the entries 

of A;, as well as, in the entries of A/+i; consequently Eq. [36] necessarily has the form 



„ /A even ^ ^ ^ ^odd , „ /rj-even ^ ^ \ odd 
Ci (A , C^Q geven ) + C2 [n , C^odd' veven ) m^^i 

but Eq. [36] must be zero for = /, for which m!^'^'^ = m^^*^ = L This implies that 
Ci = —02- By permutation invariance of the local terms we can factor out (^m^*^^ — m^*^^y 
Similarly, the homogeneity and permutation invariance of implies. 



mf^ - m?dd) (g^^^, g^^^^) ^even ^ (g^^^^ g^^^^) ^even] 

The right hand side should be zero for A/+i = /, whereby we can factor out 
(mf^^^ — mfy^^); hence the right hand side of Eq. [36] becomes 

(2fc — 1) / odd ^oddA /^even ^even\ e (n n \ ('j.7\ 
— ym^ -"^1,1 j {^2 -^1,1 j Jg (Vodd'*^evenj [oi) 

where fq (Qodd' Qeven) is a homogeneous function of order four in the entries of Qq^^ as 

well as Qeven- To evaluate fq, it suffices to let A; and A^+i be projectors of rank one where 

A/ would have only one nonzero entry on the i position on its diagonal and A/+i only one 

th 

nonzero entry on the j position on its diagonal. Further take those nonzero entries to 
be ones, giving mf^^ = mf^ = and = = 1/n. Using this choice of local terms the 
right hand side of Eq. [36] now reads 
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ETr {(kf){<,.'"|®/.)'(/.®|<,r"){<,r"|)' 

k."')k'"i®/.)(/.®kr'>(''r'i)]'} (38) 

where here the expectation value is taken with respect to random choices of local eigenvec- 
tors. Equating this and Eq. 



/. (Qodd' Qeven) = n^ETr | (|gf | ® I,) [h ® \qf^'^){q^^"\ ) (39) 



To simplify notation let us expand these vectors in the computational basis \qf^) = 
Mjiia 1^1)1^2) and \qf^^^) = "^4223 1^2)^3)- The first term on the right hand side of Eq. |38l the 



classical term, is obtained by assuming commutativity and using the projector properties, 



Tr 



Tr 



\q?){q?\^i^) (/.®kr^)(^r^ 



Tr [ui^,i2Uj^j2Vj2i,,Vk2hiUj,k2\iii2i3){jik2h\ 



Tr [(m^m) (to ^)] = Tr uv { 



[uv 



\uv 



3212 

t 



III = 



1=1 



(40) 



where ||.||p denotes the Frobenius norm and Uj are the singular values of uv. The second 
term, the quantum term, is 
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Tr 



= Tr xiv (liv)^ 



V / J2'fn2 V ' 



m2j2 



J2«2 
2 



where we used the symmetry of (^v {uv)^^ = uv (uv)^ uv {uv)^ 
Now we can calculate 



i=l 



(41) 



fi (^odd' Qeven) = \ \\uv\\p - 
giving us the desired result 



uv (uv 



(42) 



-ETr 



m 



{AU^BUY - [AQ^^BQ.Y] = d{2k- 1) fedd _ ^oddA (^even _ ^even|j43) 



from which 



X E \\uv 



uv {uv)^ 



1 — p 



(44) 



ETr {A[FBUy - ETr {AQ'^BQ^Y 
ETr {AU^BUf - ETr (AQ'^BQ^ 

d{2k-l)E(^\\uv\\p- uv{uvy 

The dependence on the covariance matrix has cancelled- a covariance matrix is one whose 
element in the i,j position is the covariance between the i*^ and j*'* eigenvalue. This shows 
that p is independent of eigenvalues of the local terms which proves the universality lemma. 

Comment: To get the numerator we used permutation invariance of A and B and local 
terms, to get the denominator we used permutation invariance of Q. 



□ 
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Comment: It is interesting that the amount of mixture of the two extremes needed to 
capture the quantum spectrum is independent of the actual types of local terms. It only 
depends on the physical parameters of the lattice. 



E. The Slider Theorem and a Summary 

In this section we make explicit use of /3— Haar properties of Q and local terms. To prove 
that there exists a < p < 1 such that the combination in Eq. [H] is convex we need to 
evaluate the expected Frobenius norms in Eq. 

2 



Lemma 3. E \\uv\ 



1/d and E 



uv (uv) 



/3^[3d(d-l)+l]+2ff(3d-l)+4 
d{l3(P+2f 



, when local terms 



have fi—Haar eigenvectors. 



Proof. It is a fact that G = uxpip-, when u is uniform on a sphere, G is a x /3— Gaussian 
matrix whose expected Frobenius norm has a distribution denoted here by xpd? (similarly 
for v). Recall that E {xD = and E {xVl = h{h + 2). 



E\\uv\\y E(x/3d2 
E||-uf||p = 



The quantum case, u'^u = 
W2 are Wishart matrices. 



iK?riip 



1 



d^ 



ttE II Gi Go II T? — o 



k=l 



d^ 



-^d ipy 



1 

d' 



, similarly v'^v 

ll<-i||p 



GIG2 



(45) 



^ = J;fi , where Wi and 



E 



uv (mv)^ 



ETr {W1W2 



ETr (W1W2) 



(46) 



hence the complexity of the problem is reduced to finding the expectation of the trace of a 
product of Wishart matrices. 



ETr {WiW2f = ETr {W1W2W1W2 



l<ijkl<d 



(47) 
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Notation 


Type 


Count 


X 


i^k & j / / 


£ {d-if 


Y 


or 

i^k k j=l 


2d^ {d - 1) 


Z 


i = k & j = I 


d^ 



Table V: Expectation values. 



where 11 denotes the product of the elements of the matrix. There are three types of 
expectations summarized in Table |Vl 
In Table M 



X 
Y 



E [n (xiXk) {yiVi)] 

2 / N 2 



E 

Z = E 

We now evaluate these expectation values. We have 

/ 

X = n 



Xi3d gfs 

Xl3{d-l) 





J 



9p 9p 
DC DC 



by QR decomposition, where gp and Xh denote an element with a /3— Gaussian and Xh 
distribution respectively; DC means "Don't Care". Consequently 



X 



n 



Xp{d-i) 





( a h\ 







n 


aXfid 913(1 + Xf}id-i)C 


= n 


axi3d 9130. 




bx^d gph + xp{d-i)d 




bxpd 9pb 
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where we denoted the four independent Gaussian entries by a,b,c,d to not confuse them as 
one number. From Eq. |17]we have 



Y = E 



4yj) (4yi 



(3d {(3d + 2) (3^ 



E {x^vY = E (xj,) E (xj) = (3d (/3d + 2) (3 {(3 + 2) . 



Eq. |l6] now reads 



E 



UV [UV] 



(3'^ [U (rf - 1) + 1] + 2(3 {3d - 1) + 4 
d {(3d^ + 2f ■ 



(48) 
□ 



Theorem. (The Slider Theorem) The quantum kurtosis lies in between the classical 
and the iso kurtoses, ^ 72 — 72- Therefore there exists a < p < 1 such that ^\ = 
PI2 + (1 ~ P) 72"*"- Further, hmAr^ooP = 1- 

Proof. We have |l - ^ES=i'?»^j} - 0' ^^^^^ T^ijlfj < Y^ijlij = ^- The last inequality 
follows from qfj < 1 . Therefore, Eq. 120) is 



7^° - 72 



a 



4 \ ^ 



- mf^^) X 



km in 



1^^) ( 

[mE (gfj - 1} < 



n (m — 1) 

From Eqs. Hi] and |40] and using the fact that the singular values cTj < 1 we have 



UV [UV] 



I l|2 

I^^IIf 



i=l 



i=l 



which proves 72 ~ 72 ^ 0- order to establish 72*° < 72 ^ 72? need to show that 



72 - 7| < 72 - 72 
HHl HE] reads 



. Eq. HH after substituting mE {qfj) = .^t^^. from Table HV] and Eqs. 



l_p = (1 _ ^2^-1) 



fc- 1 



1 - rf- 



-2fc+l 



/3 (d^ + - 2d + 1) + 4ci - 2 



(rf 



{(3d^ + 2) 
(49) 
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We want to show that < 1 — p < 1 for any integer k > 1, d > 2 and /3 > 1. All the 
factors are manifestly > 0, therefore 1 — p > 0. The first two factors are clearly < 1 so we 
need to prove that the term in the braces is too. Further, k = 1 provides an upper bound 



as 1^1 — \^ij^y2 j — i+/3d2/2 j • rewrite the term in the braces 

d {f3d^ + 2) [P + cl^ -2d+l)+Ad-2] 

{Pd^ + 2f{d+lf{d-l) ' ^ ^ 

but we can subtract the denominator from the numerator to get 



{(3d + 2) [(3 (rf^ - 2d^) + 2 (rf^ - - 1)] > V d > 2. 

This proves that flSOl) is less than one. Therefore, the term in the braces is less than one 
and hence < p < 1. Let us note the following limits of interest (recall N — 1 = 2k) 

lim (1 - p) = „ ^= 1 

1 

lim (1 - _p) ~ TT 

the first limit tells us that if we consider having two local terms and take the local dimension 
to infinity we have essentially free probability theory as expected. The second limit shows 
that in the thermodynamical limit (i.e., — i- oo) the convex combination slowly approaches 
the classical end. In the limit where /3 — )■ oo the /3 dependence in (1 — p) cancels out. This 
is a reconfirmation of the fact that in free probability theory, for /3 — )■ cxd, the result should 
be independent of {3. We see that the bounds are tight. 
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p Vs. N 




Figure 7: An example: (3 = 1: the quantum problem for all d lies in between the iso {p = 0) and 
the classical (p = 1). 

□ 

Comment: Entanglement shows itself starting at the fourth moment; further, in the 
expansion of the fourth moments only the terms that involve a pair of local terms sharing 
a site differ. Note that when the QMBS possesses a translational symmetry, there is an 
additional complication introduced by the dependence of the local terms. Though, in this 
case, the non-iid nature of the local terms complicates the matter theoretically, we have 
not seen a practical limitation of IE in our numerical experiments. 

Comment: One could from the beginning use free approximation instead of isotropic 
(m oo), in which case the proofs are simplified. 

We now summarize the main thesis of this work. We are interested in the eigenvalue 
distribution of 

^ = ^odd + ^even = Yl I ® ^M+i ® ^ + ^ ® ^M+i ® I. 

(=1,3,5, ■■■ i=2,4,6,--- 

which in a basis in that -f^odd diagonal reads H = A + Q^^BQq. Since this problem has 
little hope in being solved exactly we consider along with it two known approximations: 
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H ^ A + Q-'BQg 



Hi: 



A + Q-^BQ. 



We proved that the first three moments of the three foregoing equations are equal. We 

then calculated their fourth moments as encoded by their kurtoses (72 's) analytically and 
proved that there exists a < p <1 such that 

7|=P72' + (l-p)7r- 

It turned out that the only terms in the expansion of the fourth moments that were relevant 



were 



1 — p 



ETr |(An-ifin)' - {AQ-'BQgf^ 



(51) 



KIt{{AU-^BUy - {AQ-^BQY} 
Through direct calculation we found that the numerator ETr |(An~^Sn)^ — (^AQ^^BQgY^ 



evaluates to be 



I l|2 



uv (uv 



d {2k - 1) (mOdd - mO?d) (mf - m^J^^) E 
and the denominator ETr ^{AU-^BUf - {AQ-^BQf^ 

f^odd ^odd\ /^even ^even\ f ~ 1) V mW (n'^W 

[^2 -^1,1 j -^1,1 j I / _ ) 



1) 

Therefore 1—p does not depend on the local distribution and can generally be expressed 



as 



c^(2A; - 1)E \\uv 



1—p 



uv [UV] 



If we further assume that the local eigenvectors are /3 — Haar distributed we get 
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r/3(rf3 + ^2 _ 2rf+ 1) +4d-2" 
^ _ {d - 1) + 2) _ ■ 

Next we asserted that this p can be used to approximate the distribution 

dv'i « dv^'^ = pdv^ + (1 - p) dv''°. 

We argued that the spectra obtained using Isotropic Entanglement (IE) are accurate well 
beyond four moments. 

For illustration, we apply IE theory in full detail to a chain with Wishart matrices as 
local terms. Other types of local terms (e.g. GOE, random ±1 eigenvalues) can be treated 
similarly; therefore in Section |VT] we show the plots comparing IE with exact diagonalization 
for these cases. 



1 - d 



-2k+l 



(3dy2 



d + 1 



V. A DETAILED EXAMPLE: WISHART MATRICES AS LOCAL TERMS 

As an example take a chain with odd number of sites and for the local terms in Eq. [1] 
pick = W^W, where W is a rank r matrix whose elements are picked randomly from 
a Gaussian distribution (/3 = 1); these matrices W'^W are known as Wishart matrices. 
Clearly the maximum possible rank is r = c?^ for each of the local terms. 

Any cumulant is equal to the corresponding cumulant of one local term, denoted by 
K, times the number of summands in Eq. [H In particular, the fourth cumulant of H is 
^(N 1) _ _ 1^^^ Below we drop the superscripts when the quantity pertains to the 
whole chain. Next we recall the definitions in terms of cumulants of the mean (/x), the 
variance (o"^), the skewness (71), and the kurtosis (72) 

^ = ^2 ^ ^ ^ ^2 = ^ = ^ _ 3 . (52) 



A. Evaluation of p = ^ — 

The moments of the local terms are obtained from MOPS 



25|; 
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(53) 



mi = /3r 
m2 = f3r[f3 {r + n - 1) + 2] 
ms = /3r {(3^ [n^ + (r - 1) (3n + r - 2)] + 6/3 (n + r - 1) + 8} 
m4 = /3r {48 + [n^ + 6n'^{r - 1) + n (6r - 11) (r - 1) - 6 (r^ + 1) + + llr] 
+2/32 [6 (n2 + r^) + 17 (n(r - 1) - r) + 11] + 44/3 (n + r - 1)} 

mi^i = /3^r (r — 1) 

which for real matrices /3 = 1 yields 



rrii = r 
777,2 = r (r + + 1) 
m^ = r {n^ + 3?2 + 3rn + 3r + + 4) (54) 
1714 = r (6n^ + 21?7, + 6rn^ + 17™ + 21r + Qnr'^ + 6r^ + + + 20) 

mi_i = r (r — 1) . 

The mean, variance, skewness, and kurtosis are obtained from the foregoing relations, 
through the cumulants Eq. [52l We drop the superscripts when the quantity pertains to the 
whole chain. Therefore, using Eq. [52], we have 



fi = {N-l)r = r{N -l){n + l) 

_ n'^+Sn+i (c) _ ra^(n+6)-T-n(n+l)+21n+2r+20 



(55) 



From Eq. [22] we readily obtain 



— ETr (AU^BU) ^ = r^P (rk + n + 1?. (56) 

By The Matching Three Moments theorem we immediately have the mean, the variance 
and the skewness for the isotropic case 

;7=(A^-l)r =r{N -l){n + l) 

^ n^+3n+4 
'1 (n+l)3/2^r(Ar-l)' 

Note that the denominator in Eq. [H] becomes. 



,c ^iso '^4 '^4 



(c) _ ^(iso) 2 E {Tr {Ali^BIif - [AQ^BQY } 
72 -72 - —4 -- r2 [N-lf{n + lf 



(57) 
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In the case of Wishart matrices, m^*^*^ = mf^^^ = r, and m^^'^ = mf^^^ = r (r + n + 1), 
/T^odd _ ^even ^ ^ _ given by Eqs. [22] and [23] respectively. Therefore we can substi- 
tute these into Eq. [30] 



— E<iTr 
m 



{AlFBuY-{AQ^BQy 



TTio — m 



/3 (m — 1) / A;m (n — 1) 



(m/3 + 2) \ n (m — 1) 
Pk'^m'^ (n - 1? , 
(mp + 2) (m — 1) 



(m2 - 



(58) 



One can also calculate each of the terms separately and obtain the same results (see Appendix 
for the alternative). 
From Eq. [l3] we have 



1e 

m 



Tr {AlfBIif - Tr {AQ'^^BQ,)' 



d {2k — 1) {ni2 — fni^iY 



X E I ll^f lip — 



uv[uv] 



X 



{2k — 1) (m2 — rai l 
1 + d (d^ + d - 3) 



{<P + 2)' 

We can divide Eq. [59] by Eq. [58] to evaluate the parameter p. 



(59) 



B. Summary of Our Findings and Further Numerical Results 



We summarize the results along with numerical experiments in Tables I VI I and IVIIII to 
show the equivalence of the first three moments and the departure of the three cases in their 
fourth moment. As said above, 



Qc = V Q,so = Q Haar x Q,= [Q^^^fq^''^ 



(60) 



rp 

where, [Qf") ^ is given by Eq. [11] In addition, from Eq. [55] we can define A to be 
the part of the kurtosis that is equal among the three cases 



A = 72 



c _ 1 (rk+n+lf 
'^2 2 {n+\f 
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/3 = 1 Wishart 


Iso 


Quantum 


Classical 


Mean fi 


r{N -1) 


Variance cr^ 


r{N-l) {(f + 1) 


Skewness 71 


Vr(Af-l)(rf2+l)^ 


iE \Tr{AQlBQ.f] 


m^mf - Eq. [58] 


m^mf- Eq. [59] 


r^/c^ (r/c + n + 1)^ 


Kurtosis 72*'' 




;^E \Tt {AQ^BQg)']+A 


Eq. [55] 



Table VI: Summary of the results when the local terms are Wishart matrices. The fourth moment 
is where the three cases differ. 

Using A we can obtain the full kurtosis for the iso and quantum case, and therefore (see 
Table I VI I for a theoretical summary) : 



P 



(61) 
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Experiments based on 500000 trials 


AT Q 

iV — O 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean n 


8 


8.007 


8.007 


7.999 


Variance 


40 


40.041 


40.031 


39.976 


Skewness 71 


^VW = 1.01192 


1.009 


1.009 


1.011 


Kurtosis 72 


§1 = 0.590 


1 = 0.660 

50 


i = 0.960 


0.575 


0.645 


0.953 


Experiments based on 500000 trials 


N = 5 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean jj, 


16 


15.999 


15.999 


16.004 


Variance cr^ 


80 


79.993 


80.005 


80.066 


Skewness 71 


^^5 = 0.716 


0.715 


0.715 


0.717 


Kurtosis 72 


^ = 0-087 


^ — 0-255 


i = 0-48 


0.085 


0.255 


0.485 


Experiments based on 300000 trials 


N = 7 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean /x 


24 


23.000 


23.000 


24.095 


Variance cr^ 


120 


120.008 


120.015 


120.573 


Skewness 71 


^^30 = 0.584 


0.585 


0.585 


0.588 


Kurtosis 72 


16904 _ n nco 
206375 0.0«2 


ij = 0-153 


i = 0-320 


-0.079 


0.156 


0.331 


Experiments based on 40000 trials 


N = 9 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean // 


32 


32.027 


32.027 


31.777 


Variance a' 


160 


160.07 1 


160. 0-19 


157.480 


Skewness 71 


^VTO = 0.506 


0.505 


0.506 


0.500 


Kurtosis 72 


539142 _ Q^g4 
3283175 u.j.u'i 


= 109 

800 


i = 0-240 


-0.165 


0.109 


0.213 


Expcnimciits l)as(Hi 011 2000 (rials 


N = 11 


Theoretical value 


Numerical Experiment 




Iso 


Quantum 


Classical 


Iso 


Quantum 


Classical 


Mean /j, 


40 


39.973 


39.973 


39.974 


Variance 


200 


200.822 


200.876 


197.350 


Skewness 71 


^V2 = 0.452548 


0.4618 


0.4538 


0.407 


Kurtosis 72 


11162424 n 01 q 
52454375 U.zio 


m = 0-084 


m = 0-192 


-0.189 


0.093 


0.102 



Tabic VIII: The mean, variance and skewness of classical, iso and quantum results match. However, 
the fourth moments (kurtoses) differ. Here we are showing results for d = 2, r = 4 with an accuracy 
of three decimal points. 37 



The numerical convergence of the kurtoses to the theoretical values were rather slow. To 
make sure the results are consistent we did a large run with 500 million trials for = 5, 
d = 2, r = 3 and /3 = 1 and obtained four digits of accuracy 



72 - Yr 



0.39340 Numerical experiment 
0.39347 Theoretical value. 
Convergence is faster if one calculates p based on the departing terms alone (Eq. IFIi) . In 



Y2 - ir 



this case, for full rank Wishart matrices with N = 5 and d = 2 



/3 = 1, trials: 5 Million 


1 — p 


Numerical Experiment 


0.57189 


Theoretical Value 


0.57183 



/3 = 2, trials: 10 Million 


1 — p 


Numerical Experiment 


0.63912 


Theoretical Value 


0.63938 



Below we compare our theory against exact diagonalization for various number of sites 
A^, local ranks r, and site dimensionality d (Figures [8l fT2ll . 
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N=3, d=2, r=4, trials=500000 
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• Exact Diagonalization 
— I.E. Theory: p = 0.07 
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Figure 8: = 3 examples. Note that the last two plots have the same p despite having different 
ranks r. This is a consequence of the Universality Lemma since they have the same and d. 
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N=5, d=2, r=4, trials=500000 
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N=5, d=3, r=5, trials=2000 
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Figure 9: = 5 
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— I.E. Theory: p = 0.59 
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Figure 10: iV = 7 
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N=9, d=2, r=2, trials=40000 
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— I.E. Theory: p = 0.68 
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Figure 11: = 9. Note that the two plots have the same p despite having different local ranks. 



N=11,d=2, r=4, trials=2800 
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— I.E. Theory: p = 0.73 




Figure 12: iV = 11 
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VI. OTHER EXAMPLES OF LOCAL TERMS 



Because of the Universality lemma, p is independent of the type of local distribution. 
Furthermore, as discussed above, the apphcation of the theory for other types of local terms 
is entirely similar to the Wishart case. Therefore, we only show the results in this section. 
As a second example consider GOE's as local terms, i.e., = ^-j^, where G is a full 

rank matrix whose elements are real Gaussian random numbers. 
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Figure 13: GOE's as local terms 



Lastly take the local terms to have Haar eigenvectors but with random eigenvalues ±1, 
i.e., -f^z,/+i = QfAiQi, where A; is a diagonal matrix whose elements are binary random 
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N=3, d=5, trials=50000 
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Figure 14: Local terms have a random binomial distribution. 

variables ±1 (Figure [H]). 

In this case the classical treatment of the local terms leads to a binomial distribution. As 
expected p = 1 in Figure [H] has three atoms at —2, 0, 2 corresponding to the randomized 
sum of the eigenvalues from the two local terms. The exact diagonalization, however, shows 
that the quantum chain has a much richer structure closer to iso; i.e, p = 0. This is captured 
quite well by IE with -p = 0.046. 



VII. BEYOND NEAREST NEIGHBORS INTERACTION: L>2 

If one fixes all the parameters in the problem and compares L > 2 with nearest neighbor 
interactions, then one expects the former to act more isotropic as the number of random 
parameters in Eq. [1] are more. When the number of random parameters introduced by the 
local terms, i.e., (A^ — L + 1) and c?^ are comparable, we find that we can approximate 



the spectrum with a high accuracy by taking the summands to be all isotropic 29|] (See 
Figures [IMIZ!). 
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Figure 15: IE method approximates the quantum spectrum by H^^ = X]f=i Qf Hi^... ^i^2Qi 
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Figure 16: IE method approximates the quantum spectrum by H^^ = Ylf=i Qf ^i,- - ,i+3Qi- 
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Figure 17: IE method approximates the quantum spectrum by H 



Td=iQlHi,-,i+iQl 



Most distributions built solely from the first four moments, would give smooth curves. 
Roughly speaking, the mean indicates the center of the distribution, variance its width, 
skewness its bending away from the center and kurtosis how tall and skinny versus how short 
and fat the distribution is. Therefore, it is hard to imagine that the kinks, cusps and local 
extrema of the quantum problem (as seen in some of our examples and in particular Figure 
in [T4|) could be captured by fitting only the first four moments of the QMBS Hamiltonian to 
a known distribution. It is remarkable that a one parameter (i.e., p) interpolation between 
the isotropic and classical suffices in capturing the richness of the spectra of QMBS. 

VIII. CONJECTURES AND OPEN PROBLEMS 

In this paper we have offered a method that successfully captures the density of states of 
QMBS with generic local interaction with an accuracy higher than one expects solely from 
the first four moments. We would like to direct the reader's attention to open problems that 
we believe are within reach. 

1. We conjecture that the higher moments may be analyzed for their significance. For 
example, one can show that the fraction of departing terms in the expansion of the 
higher moments (e.g. analogous to bold faced and underlined terms in Eqs. [T6IT8] 
but for higher moments) is asymptotically upper bounded by In Sectior UTCl we 

conjectured that their expectation values would not change the moments significantly. 
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It would be of interest to know if 

ett{...q-^b-^qa-^q-^b^^q...} < 

ETr {. . . Q-'B^'Q,A^'Q-'B^'Q, . . .} < 
ETr{...Il-^B-^IlA^^U'^B-^U...} . 
For example, one wonders if 

ETr I [AQ-^BQY^ < ETr | [AQg^BQ^Y^ < ETr | [AU-^BuY^ 

for > 2; we have proved that the inequality becomes an equality for k = 1 (Departure 
Theorem) and holds for k = 2 (Slider Theorem). 

2. Though we focus on decomposition for spin chains, we believe that the main theorems 
may generalize to higher dimensional graphs. Further rigorous and numerical work in 
higher dimensions would be of interest. 

3. At the end of Section HVDI we proposed that more general local terms might be treated 
by explicitly including the extra terms (Type I terms). 

4. Application of this method to slightly disordered systems would be of interest in con- 
densed matter physics. In this case, the assumption of fully random local terms needs 
to be relaxed. 

5. In our numerical work, we see that the method gives accurate answers in the presence 
of an external field. It would be nice to formally extend the results and calculate 
thermodynamical quantities. 

6. We derived our results for general (3 but numerically tested /3 = 1, 2. We acknowledge 
that general (3 remains an abstraction. 

7. Readers may wonder whether it is better to consider "iso" or the infinite limit which 
is "free". We have not fully investigated these choices, and it is indeed possible that 
one or the other is better suited for various purposes. 

8. A grander goal would be to apply the ideas of this paper to very general sums of 
matrices. 
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IX. APPENDIX 



To help the reader with the random quantities that appear in this paper, we provide 
explanations of the exact nature of the random variables that are being averaged. A common 
assumption is that we either assume a uniformly randomly chosen eigenvalue from a random 
matrix or we assume a collection of eigenvalues that may be randomly ordered, Random 
ordering can be imposed or a direct result of the eigenvector matrix having the right property. 
Calculating each of the terms separately and then subtracting gives the same results. 



Ik 

m 



Tr{AQ^BQy 




^ ] Cli J {^hjl^hj2^i2jl^i2j2) \ 5 (62) 

where a, and bj are elements of matrices A and5 respectively. The right hand side of Eq. 
[62] can have terms with two collisions (i.e., ii = 12 and ji = 72), one collision (i.e. ii 7^ 22 
exclusive-or ji 7^ J2), or no collisions (i.e., ii 7^ 12 and ji 7^ j2)- Our goal now is to group 
terms based on the number of collisions. The pre-factors for two, one and no collisions along 
with the counts are summarized in Table IIV[ Using the latter we can sum the three types 
of contributions, to get the expectation 



Tr(AQ^fig)' 



(/3+2) 
(m/3+2) 



E (a2) E (62) + 



^ [E (62) E (aias) + E (a') E (6162)] - £tSe ia,a2) E (6162) • 



(63) 



(m^+2) V" ^ V"-^""^/ ' \^ J \"l^2Ji (m/3+2)' 

If we take the local terms to be from the same distribution we can further express the 
foregoing equation 



— E 
m 



Tr (AQ^Sg)' 



1 



(m/3 + 2] 



[iP + 2) m2 + /3 (m - 1) E (0102) {2m2 - E (0102)}] . 



(64) 



The quantity of interest is the difference of the classical and the iso (see Eq. [2U]) . 



{mp+2) ^ ' 
(m/3+2) 



iETr {ABf - ^ETr {AQ^ BQ)' = 
E (62) - E (62) E (0102) - E (a2) E (6162) + E (0102) E (6162)} = 

' (A) (A) (B) . (A) ^R^^ 



, (A) (B) (B) {A 

I 7712 "^2 ~ ''^2 "^1,1 

(^{A) _ (A) 
(m/3+2) y'h ''h,l 



-771)2 
( (' 



(B) {B)\ 

" — 7n\^^ 1 = 



^1,1 / 



(65) 
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If we assume that the local terms have the same distribution: m2 = m!'^^ = rn"^\ 
mil = ^^11 = rn[^^ as in Eq. [6ll the foregoing equation simplifies to 



lETr {ABf - -ETt [AQ^BQ] ' = ^^^^ {m, - mi,if . 
m m ^ ' [mp + 2) 

In the example of Wishart matrices as local terms we have 



mi^i = E (aia2) = k (k 



1712 = E (a^) = rk {rk + n + 1) 
1) r2 + ^ { (tn^-i - 1) (n + r + 1) + tn''-^ {n - 1) (r - 1) } 



(fc - 1) + {tn''~^{nr + 2)-n-r-l}. 



(66) 
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